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Abstract 

Since the end of 80's, the region of sky of galactic coordinates (Z ~ 309°, b ~ 18°), corre- 
sponding to a declination 7 ~ —44° and right ascension a ~ 202°, usually denoted as the 
"Great Attractor", is known to control the overall galaxy flow in our local Universe. In 
this sense, this direction might represent a natural candidate to characterize a hypothetical 
Earth's "absolute motion". Our analysis of the extensive ether-drift observations recently 
reported by an experimental group in Berlin provides values of a and 7 that coincide almost 
exactly with those of the Great Attractor and not with the values 7 ~ —6° and a ~ 168° 
obtained from a dipole fit to the anisotropy of the CMB. This supports in a new fashion the 
existence of a discrepancy between the observed motion of the Local Group and the direction 
obtained from the CMB dipole. 
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1. Introduction 



Since the end of 80's, a region of the sky in the direction of the Centaurus cluster, the so called 
"Great Attractor", at galactic coordinates (7 ~ 309°, b ~ 18°), corresponding to a declination 
j ~ —44° and right ascension a ~ 202°, is known to play an important role in describing the 
deviations of galaxies from a pure Hubble flow in our local Universe, ^ HJ HI 03 E] • In this 
sense, the Great Attractor, marking a preferred direction in space, might represent a natural 
candidate to characterize a hypothetical Earth's "absolute motion". 

On the other hand, since the discovery of an anisotropy in the cosmic microwave back- 
ground (CMB), it has been generally accepted that the kinematical parameters for such an 
absolute motion should coincide with those deduced from a dipole fit to the COBE data 
j7j. In this context, one predicts a velocity of the Solar System v ~ 370 km/s, an average 
declination angle 7 ~ —6° and a right ascension a ~ 168°. These values have usually been 
adopted in the interpretation of the data from the ether-drift experiments in a laboratory. 

However, the latest ether-drift experiments, combining the possibility of active rotations 
of the apparatus with the use of cryogenic optical resonators El EH ; have reached such a 
high precision (0(1O~ 16 ) in the relative frequency shifts) to require a fully model-independent 
analysis of the data. In fact, the physical nature of a hypothetical preferred frame is still 
unknown. Therefore, assuming from the very beginning one particular set of values for 
(w,7,a) one might introduce uncontrolled errors in the interpretation of the experimental 
results. 

This is even more true noticing that a fully model-independent analysis of the ex- 
tensive ether-drift observations reported by Herrmann et al. in Ref.0 provides an average 
(absolute) value of the declination angle I7I ~ 43° ± 3° that would rather favour an alterna- 
tive of the type represented by the Great Attractor. In this paper we'll further extend the 
analysis of Ref. jllj that, being limited to a restricted set of observables, could not determine 
the value of a and the sign of 7. As we shall illustrate, our new results, after inclusion of 
other observable quantities from Ref.t9;, provide values of a and 7 that are in remarkable 
agreement with those of the Great Attractor. For this reason, our results, while providing the 
first modern evidence for an ether drift from a laboratory experiment, support the indications 
obtained from the observed motion of galaxies. 

The plane of the paper is as follows. In Sect. 2 we shall report the relevant formalism 
used in the analysis of the ether-drift experiments and the basic experimental data of Ref. [3]. 
In Sect. 3 we shall present our analysis of these data while in Sect. 4 we shall summarize our 
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results and present our conclusions. 



2. General formalism and experimental data 

The starting point for our analysis is the expression for the relative frequency shift of the two 
optical resonators at a given time t. This is expressed as (HI 

—Q- = S{t) sin 2uj Iot t + C{t) cos 2uj Iot t (1) 

where uj rot is the rotation frequency of one resonator with respect to the other which is 
kept fixed in the laboratory and oriented north-south. The Fourier expansions of the two 
amplitudes S(t) and C(t) are predicted to be 

S(t) = S s i sin r + S c i cos r + S S 2 sin(2r) + S C 2 cos(2r) (2) 

C(t) = C + C s i sin r + C c i cos r + C s2 sin(2r) + C c2 cos(2r) (3) 

where r = w s id^ is the sidereal time of the observation in degrees and w s id ~ 23^56' ' 

Introducing the colatitude of the laboratory x(~ 37.5° for Berlin), one finds the expres- 
sions reported in Table I of Ref. [S], 

Co = -i^^(3cos2 7 -l), (4) 
o 

C s i = -if sin 27 sin a sin 2x, (5) 

C c i = -if sin 2 7 cos a sin 2x, (6) 

C S 2 = ^-if cos 2 7 sin 2a (1 + cos 2 x), (7) 

C c2 = iifcos 2 7 cos2a(l + cos 2 x) (8) 
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where 

K = (1/2- + 5)^ (9) 

cr 

and (1/2— (3+6) indicates the Robertson-Mansouri-Sexl ^1^3 (RMS) anisotropy parameter. 
The corresponding S— quantities are also given by S s \ = — C c \/ cosx, S c i = C s i/cosx> 
S s2 = — -| 2cos ^ C c2 and S C 2 = 2cos x Q „ it might be interesting that these relations can 

*•* l+cos^ x 1+cos^ x 

also be derived from an old paper by Nassau and Morse, published in the Astrophysical 
Journal about eighty years ago (see Eqs. (20-24) of Ref.|l4j. 



2 



Table 1: The experimental values of the C— coefficients for the 15 observation periods of 
Ref.gj. 



Observation i C s i[xl0~ 16 ] C c i[xl0~ 16 ] C^fxlO" 16 ] C c2 [xl0- 16 ] 



1 


-2.7 ±4.5 


5.3 ±4.8 


-3.2 ±4.7 


1.2 ±4.2 


2 


-18.6 ±6.5 


8.9 ±6.4 


-11.4 ±6.5 


-5.0 ±6.4 


3 


-0.7 ±3.9 


5.3 ±3.6 


5.0 ±3.5 


1.6 ±3.8 


4 


6.1 ±4.6 


0.0 ±4.8 


-8.1 ±4.8 


-4.0 ±4.6 


5 


2.0 ±8.6 


1.3 ±7.7 


16.1 ±8.0 


-3.3 ±7.2 


6 


3.0 ± 5.8 


4.6 ± 5.9 


8.6 ± 5.9 


-6.9 ± 5.9 


7 


0.0 ± 5.4 


-9.5 ± 5.7 


-5.5 ±5.6 


-3.5 ±5.4 


8 


-1.1 ±8.1 


11.0 ±7.9 


0.9 ±8.3 


18.6 ±7.9 


9 


8.6 ±6.5 


2.7 ±6.7 


4.3 ±6.5 


-12.4 ±6.4 


10 


-4.8 ±4.8 


-5.1 ±4.8 


3.8 ±4.7 


-5.2 ±4.7 


11 


5.7 ±3.2 


3.0 ±3.4 


-6.3 ±3.2 


0.0 ±3.5 


12 


4.8 ±8.0 


0.0 ±7.0 


0.0 ±7.6 


1.5 ±7.7 


13 


3.0 ±4.3 


-5.9 ±4.3 


-2.1 ±4.4 


14.1 ±4.3 


14 


-4.5 ±4.4 


-2.3 ±4.5 


4.1 ±4.3 


3.2 ±4.3 


15 


0.0 ±3.6 


4.6 ±3.4 


0.6 ±3.2 


4.9 ±3.3 



Notice a critical detail for the separation of the signal in its elementary components. A 
small mismatch in the definition of the sidereal time, say r true = r + At, induces a rotation 
of the various parameter pairs of angles At and 2Ar with a corresponding re-definition of 
the right ascension a true = a + At. 

The experimental data reported in Ref.jH] were obtained during 15 short-period observa- 
tions performed from December 2004 to April 2005. As suggested by the same authors, it 
is safer to concentrate on the observed time modulation of the signal, i.e. on the quantities 
C s i, C c i, C S 2, C C 2 and on their 5-counterparts. In fact, the constant components C = C$ and 
S = So are most likely affected by spurious systematic effects such as thermal drift (see also 
the discussion in Ref. [S]). 

The individual determinations of the various parameters, for each of the 15 short-period 
observations, as extracted from Fig. 3 of Ref.[S], are reported in Table 1 and Table 2. 
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Table 2: The experimental values of the S— coefficients for the 15 observation periods of 
Ref.0. 



Observation i 


Ssi [xl(T 16 ] 


S c i[xl0- 16 ] 


S s2 [xl0- 16 ] 


S c2 [xl0- 16 ] 


1 


11.2 ±4.7 


11.9 ±4.9 


1.8 ±4.9 


0.8 ±4.5 


2 


1.8 ±6.5 


-4.3 ±6.5 


6.4 ±6.4 


1.8 ±6.4 


3 


-3.3 ±3.8 


2.9 ±3.8 


-5.9 ±3.8 


4.6 ±4.0 


4 


12.7 ±5.1 


14.3 ±5.5 


-1.9 ±5.3 


-3.3 ±5.1 


5 


4.7 ±8.4 


-6.9 ±7.3 


-1.8 ±8.0 


-7.8 ±7.0 


6 


5.2 ±5.8 


-3.0 ±5.9 


7.1 ±5.9 


-5.9 ±5.8 


7 


11.1 ±5.3 


-13.4 ±5.4 


-4.5 ±5.5 


-9.8 ±5.5 


8 


—12.1 ± 8.9 


0.0 ± 8.8 


-3.1 ± 9.0 


1.4 ± 8.9 


9 


-4.8 ±6.3 


6.5 ±6.4 


-8.1 ±6.3 


3.5 ±6.5 


10 


9.8 ±5.0 


4.8 ± 5.0 


1.9 ±5.0 


-9.2 ±4.8 


11 


0.0 ±3.2 


-3.9 ±3.6 


1.0 ±3.1 


-2.2 ±3.4 


12 


-12.7 ±7.7 


8.5 ±6.8 


-8.3 ± 7.2 


-7.1 ±7.4 


13 


-7.9 ±4.7 


-4.3 ±4.8 


-1.9 ±4.8 


-6.2 ±4.7 


14 


16.1 ±4.9 


12.0 ±5.2 


2.9 ±4.9 


-9.6 ±4.8 


15 


13.9 ±3.9 


-7.0 ±3.4 


-3.3 ±3.5 


3.0 ±3.6 



3. Analysis of the data 

The analysis of Ref . was restricted to the combinations 



Cll = A 




(10) 


C22 = A 


/Cs2 + Cc2 


(11) 


Sn = 1 


J^sl + $cl 


(12) 


•522 = 1 


J Sg 2 + 5*c2 


(13) 



This was useful to reduce the model dependence in the analysis of the data. In this way, 
in fact, the right ascension a and any possible uncertainty related to the definition of the 
sidereal time drop out from the theoretical predictions that will only depend on (7!, and the 
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Table 3: The experimental values of the combinations of C— and S— coefficients defined in 
Eas.(|lU j) -(|13 jl as obtained from our Table 1 and Table 2. 



VJUdCI V cLXjLKJLl I 


UH|X1U J 


r'nahci n _16 i 




QonTvi n _16 i 

022|XxU J 


1 

± 


O.a Zl <±. . 


^ ^ 4 4 fi 

O.O ZL 4:.U 


±U.O ZL *±.o 


9 D 4- 4 Q 


o 
z 


iU.U HZ U.4: 


i o K _i_ a z. 

li.t) ZL U.U 


A ft _|_ ft K 

4:.U ZL U.U 


U.U ZL u.*± 


Q 
O 


r; q _i_ q a 
O.O ZL o.u 


O.O ZL o.u 


^.^A- ZL O.O 


1 .(J ZL O.O 




« 1 _|_ 4 R 

U. X ZL 4:. U 


C/.U ZL ^±.0 


i q i _l c o 


O.O ZZ . X 


'J 


.1 ZL O . 4: 


1ft K _|_ q n 

xu.u zz o.u 


0.4: ZL / . / 


O.u ZL 1 . X 


u 


u. u zc o.tj 


11 n 4- ^ Q 

X J. .U Z[Z u. <J 


ft n _i_ r: q 

U.U ZL (J.C? 


q o _l k q 

U.L* ZL (J.C/ 


7 


9.5 zt 5.7 


6.5 ±5.5 


17.4 ± 5.4 


10.7 ±5.5 


8 


11.0 ±7.9 


18.7 ± 7.9 


12.1 ±8.9 


3.4 ±9.0 


9 


9.1 zb 6.5 


13.1 ±6.4 


8.1 ±6.4 


8.8 ±6.4 


10 


7.0 ±4.8 


6.5 ±4.7 


10.9 ± 5.0 


9.4 ±4.8 


11 


6.4 ±3.1 


6.3 ±3.2 


3.9 ±3.6 


2.4 ±3.4 


12 


4.8 ±8.0 


1.5 ±7.7 


15.3 ± 7.4 


10.9 ±7.3 


13 


6.6 ±4.3 


14.3 ±4.3 


9.0 ±4.7 


6.5 ±4.7 


14 


5.1 ±4.5 


5.2 ±4.3 


20.0 ± 5.0 


10.0 ±4.8 


15 


4.6 ±3.4 


5.0 ±3.3 


15.6 ±3.8 


4.4 ±3.5 



overall normalization \K\. The relevant numbers for these auxiliary quantities can be found 
in Table 3. Thus, we obtain the relations 

Cn = i|A-|sin2| 7 |sin2x C 22 = ±\K\ cos 2 7 (1 4 cos 2 X ). (14) 

The corresponding ^-coefficients are also predicted as Sn = C\\j cosx and S22 = i+c°s*x 
C22- Using the weighted averages of the values in Table 3 for these coefficients 

(C n ) = (6.7 ±1.2) ■ 10~ 16 (C22) = (7.6 ± 1.2) • 10~ 16 (15) 

(Sil) = (11.0 ± 1.3) • 10~ 16 (S 22 ) = (6.3 ± 1.3) • 10~ 16 (16) 
one gets |11| an average declination 

| 7 | ~ 43° ± 3° (17) 
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Table 4: The absolute value of the declination angle | and the normalization factor \K 
obtained from the data in Table 3 for groups of three observation periods. 

Observations I7I [degrees] |-ftT|[xlO~ 16 ] 

(1-3) 46+£ 29 ± 6 

(4-6) 32+^ 33 ± 7 

(7-9) 39tg 43 ± 7 

(10-12) 42± 7 9 27 ± 5 

(13-15) 4ll? 36 ± 5 



and an average normalization 

|i^T| ~ (33 ± 3) • (18) 

To check the regularity of the data, one can also extract the average declination angle 
from various groups of observation periods. Packing the data in groups of three observations, 
one gets the results shown in Table 4. Notice the remarkable consistency among the various 
sets of data. 

To extend the above analysis and get information on a and the sign of 7, we shall now 
try to re-construct the full amplitudes S(t) and C(t) of Eq.Q or, more precisely, their 
variable parts S(t) — So and C(t) — Co- In this way, in fact, possible problems related to the 
deconvolution of the signal in its elementary components drop out. To this end, we shall use 
Eqs.(J2J) and ©to generate a suitable signal where: 

i) the individual C— and S— coefficients are fixed to their experimental values of 
Ref.pj reported in our Tables 1 and 2 

ii) the average times for the 15 individual observation periods are taken from Fig. 3 of 
Ref.jHj. Once the figure is reproduced on the screen, these time coordinates can be extracted 
to a good accuracy from the number of pixels Npixel of the data points. In this way we get 
(in days since 1/1/2000) t = ( n p™'- 475 ) + 1900 where Npixel= 317, 344, 351, 365, 377, 400, 
412, 423, 427, 437, 470, 493, 496, 503, 507 for the 15 observation periods. 

The resulting data sets are reported in Table 5. 
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Table 5: The re-constructed amplitudes obtained following the steps i) and ii) described in 
the text. 



Observation i (C(U) - C )[xl0- 16 ] (S(U) - 5 )[xl0- 16 ] 



1 


2.7 ±6.8 


17.0 ±6.9 


2 


19.0 ±9.2 


-7.6 ±9.2 


3 


9.3 ±5.5 


0.9 ±5.5 


4 


13.3 ±6.8 


12.1 ± 7.5 


5 


-13.6 ± 12.2 


11.7 ± 11.3 


6 


-7.7 ±8.3 


-13.3 ±8.3 


7 


-7.6 ±8.1 


-20.7 ±7.8 


8 


-2.4 ± 11.7 


-3.0 ± 12.6 


9 


3.1 ±9.5 


3.4 ±9.1 


10 


12.0 ±6.8 


-6.3 ±7.1 


11 


0.6 ±4.9 


-4.5 ±4.9 


12 


-1.9 ± 11.3 


-11.1 ± 10.5 


13 


4.1 ±6.2 


-0.2 ±6.8 


14 


2.9 ±6.4 


-18.8 ±7.0 


15 


-1.4 ±5.0 


-14.7 ±5.1 



We can now try to fit the values of Table 5 to the theoretical predictions (/t = ±1) 

C(t) C - ^ K ^ s[n2 X s[n2 "f cos (® ~ t) | (l±cos 2 x)cos 2 7cos2(a-r) ^ ^ 

S(t) -S = M|K|( SinXSin27 2 Sill(a ~ T) + COSXCOs27 2 Sin2(a " T) ) (20) 

to obtain information on \K\, a and the sign of 7. 

Before presenting the results, we observe that the structure of the problem is such that 
one should first suitably constrain the fit. In fact, all together there are 16 possible choices 
arising from the 2 possible signs of fi, the 2 possible signs of 7 and the 4 possible choices of 
a (sin a = ±|sina| and cos a = ±|cosa|). Considering all possibilities we have found that 
the various fits to the values of C(t) — Cq do not provide any definite information. In fact, 
in all cases the quality of the chi-square is very close to that of the " null result" defined by 
I if I = and the parameters cannot be constrained in any meaningful way. 
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Considering the values of S(t) — So the situation is different. In fact the fit routine 
(MINUIT) finds in this case two configurations whose chi-square is definitely lower than the 
chi-square of the " null-result" (~ 38) . The existence of multiple solutions had to be expected 
since the C— (and S— ) coefficients of Eqs.(4)-(8) are unchanged under the simultaneous 
replacements 7 — > —7 and a — > a + ir. The absolute minimum, whose chi-square is ~ 16 (for 
12 degrees of freedom), corresponds to /i = — 1 with the physical parameters being in the 
following ranges 

\K\ = (39 ± 15) • 1(T 16 7 = -30°;22° " = 204° ±12° (21) 

As one can see, the values for \K\ and | 7 | from our reconstructed amplitude for Sit) — Sq are 
in good agreement with Eas.(|17|) and (|18|) . and with those reported in our Table 4, that were 
deduced from the average values of the positive-definite coefficients Cn, C22, Su and S^- 

We have no explanation for the different behaviour obtained using our re-constructions of 
C(t) — Co and S(t) — So. Perhaps this might be due to the asymmetric experimental set-up 
where only one cavity is rotated or, since the 15 data sets, each spanning from 24 hours to 100 
hours in length, have been summarized into just a single point, the inherent inexactness in 
our re-construction affects differently the two sets of data. On the other hand, the indications 
we have obtained from 5(t) — So are so consistent with the previous results of Ref . [TT] . and 
with those reported in our Table 4, to suggest that they should still persist in a more refined 
analysis using the full raw data or after inclusion of new observations. 

4. Summary and conclusions 

The possibility of a large concentration of matter in the region of sky 7 ~ —44° and a ~ 202°, 
the "Great Attractor" , was proposed to describe the deviations of galaxies from a pure Hubble 
flow [Q-[H in our local Universe. Its inclusion in multi-attractor models, see e.g. Ref. |15j. 
provides a successful representation of the observed velocity field for large samples of galaxies, 
such as the ~ 3400 ones contained in the MARK III catalog |18j . 

The resulting motion of the Local Group is known to exhibit some discrepancy with 
respect to the direction obtained from a dipole fit to the anisotropy of the CMB. The value 
of the misalignment angle is both sample- and model dependent. For instance in the multi- 
attractor model of Ref.|15j it ranges from 11° to 49° with a typical value of ~ 30°. 

In principle, the existence of this discrepancy might also show up in ether-drift experi- 
ments. Up to now, these have been analyzed assuming that the kinematical parameters of a 
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hypothetical Earth's "absolute motion" should coincide with those extracted from a dipole fit 
to the COBE data. Rather, in our opinion, one should leave out the angular variables 7 and 

2 

a and the quantity K = (1/2 — (3 + 8)\ as free parameters and check the overall agreement 
with the data. 

In this paper, using the extensive ether-drift observations reported by Herrmann et al. [§], 
we have obtained evidence for an angular pair that sizeably differs from the values a ~ 168° 
and 7 ~ —6° obtained from a dipole fit to the anisotropy of the CMB. The motion we have 
detected points toward 7 = — 2>Q°^^ and a = 204° ± 12° and the central values coincide 
with those of the Shapley supercluster [Ej. However, the error in 7 is large. Therefore, we 
could use the result of Ref.|llj I7I = 43° =L 3°, to further sharpen our estimate. In this way, 
the angular parameters of the Earth's motion coincide almost exactly with those of the core 
of the Great Attractor (7 ~ -44° and a ~ 202°). 

Actually, the agreement with the position of the Great Attractor is even too good. In 
fact, an ether-drift observation on the Earth contains, in principle, the effect of the solar 
motion ('sm') relatively to the centroid of the Local Group. For instance, using for this 
motion the values of Ref.jTTl, v sm ~ (306 ± 18) km/s, 7 sm ~ (51° ± 6°), a sm ~ (332° ±11°), 
and using for the peculiar motion of the Local Group the values of Ref.|15j obtained fitting 
the multi- attractor parameters to the full MARK III catalog, namely Vlg ~ 488 ± 80 km/s 
(in the CMB frame), 7lg ~ —50° ± 8° and olg ~ 202° ± 10°, one can predict a theoretical 
velocity through 

v th = V LG + v sm (22) 

obtaining the following values 

v th = (276 ± 71) km/s 7th = -30°t^ a th = 240°^ (23) 

These might be compared with the results of our fit (|21jl and with our Eas. (|17jl and (|18jl 
having a theoretical estimate of |(l/2 — (3 + 5)\ to transform into a velocity value, through 
Eq.@, the precise indication (J18|) on the normalization factor \K\. 

To this end, we observe preliminarily that our model-independent analysis leads to rather 
large values of the RMS anisotropy parameter. In fact, using our result \K\ ~ (33 ±3) • 10~ 16 
and the theoretical range u t h ~ (276 ± 71) km/s reported above, one obtains a range of the 
RMS parameter 25 • 10~ 10 < [(1/2 - (3 + 5)\ < 70 • 10~ 10 (with a central value 39 • 10 -10 ) in 
good agreement with the theoretical prediction | ( 1/2 — (3 + J)[ t h ~ 42 • 10~ 10 of Refs. [T§ll2T)] . 
Equivalently, using the two theoretical inputs |(l/2 — f3+5)\ t h ~ 42-10 -10 and v t h ~ (276±71) 
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km/s one predicts 20 ■ 10~ 16 < \K\ th < 56 • 10" 16 (with a central value 36 • 10" 16 ) m good 
agreement with our Eq. (|18|) , 

This should be compared with the result of Ref.|5| (where the values v ~ 370 km/s, 
a ~ 168° and 7 ~ —6° were assumed) |(l/2 — (3 + 5)\ ~ (2 ±2) ■ 10~ 10 . In this way, one would 
predict \K\ ~ (3 ±3) • 10" 16 which is one order of magnitude smaller than the value reported 
in our Ea. ()18|) . In this sense, our results show that the very small RMS parameter reported 
in Ref. 9 , rather than being due to the smallness of the signal, might originate from more 
or less accidental cancellations among the various entries. 

Acknowledgements 

We thank C. Marinoni for kindly providing the estimate of the errors on the velocity 
parameters reported in Table 6 of Ref. |15j 



10 



References 

[1] D. Lynden-Bell et al.,ApJ. 326, 19(1988). 
[2] A. Dressier, ApJ. 329, 519 (1988). 

[3] S. M. Faber and D. Burstein, in Large-Scale motions in the Universe, V. C. Rubin and 
G. Coyne Eds., Princeton University Press 1989, p. 115. 

[4] D. Burstein, Rep. Progr. Phys. 53, 421 (1990). 

[5] J. R. Mould et al., ApJ.545, 547 (2000) [Erratum]. 

[6] R. Kaldare, M. Colless, S. Raychaudhury and B. A. Peterson, Mon. Not. R. Astron. Soc. 
339, 652 (2003). 

[7] G. F. Smoot et al., ApJ.371, LI (1991). 

[8] P. Antonini, M. Okhapkin, E. Goklu and S. Schiller, Phys. Rev. A71, 050101(R) (2005) 
arXiv:gr-qc/0504109 . 

[9] S. Herrmann, A. Senger, E. Kovalchuk, H. Miiller and A. Peters, Phys. Rev. Lett. 95, 
150401 (2005) larXiv:physic s /05080971 . 

[10] P. L. Stanwix et al., Phys. Rev. Lett. 95, 040404 (2005) |arXiv:hep-p h/0506074 . 

[11] M. Consoli and E. Costanzo, arXiv:gr-qc/0511160 submitted to Phys. Rev. Letters. 



[12] H. P. Robertson, Rev. Mod. Phys. 21, 378 (1949). 

[13] R. M. Mansouri and R. U. Sexl, Gen. Rel. Grav. 8, 497 (1977). 

[14] J. J. Nassau and P. M. Morse, ApJ. 65, 73 (1927). 

[15] C. Marinoni, P. Monaco, G. Giuricin and B. Constantini, ApJ. 505 (1998) 484. 

[16] D. Proust et al., |arXiv:astro-ph/050 9903 . 

[17] S. Courteau and S. Van Den Bergh, AJ. 118 (1999) 337. 

[18] J. A. Willick et al., ApJS 109, 333 (1997). 

[19] M. Consoli, A. Pagano and L. Pappalardo, Phys. Lett. A318, 292 (2003). 

[20] M. Consoli and E. Costanzo, Phys. Lett. A333, 355 (2004). 



11 



